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LUNISOLAR PERTURBATIONS OF THE MOTION 
OF ARTIFICIAL SATELLITES 

SUMMARY 


An expression for the Juni solar disturbing function is given which 
provides the basis for the calculation of perturbations in the motion of 
artificial satellites. Previous work in this field is described in references 
1, 2, 3, and 4. The Lagrange planetary equations are derived in a form 
suitable for satellites of arbitrary eccentricity in contrast to the standard 
form which depends upon first developing the disturbing function in powers 
of eccentricity. 

The spectrum of the disturbing function is partitioned into secular, 
long period, intermediate, and short period terms. The motion of the short 
period terms is dominated by the mean anomaly of the satellite. The re- 
maining terms are derived by averaging the disturbing function over the 
satellite's mean anomaly. Those terms which contain the mean longitude 
of the sun or the moon are called terms of intermediate period. For the 
sun these terms have periods of multiples of 180 days, for the moon the 
corresponding periods are multiples of 14 days. Additional averaging over 
the mean anomaly of the disturbing body results in the secular and long 
period terms. 

The perturbations arising from the secular and long period terms 
have been discussed in references 5 and 6. In this paper analytic formulas 
forlunisolar perturbations of intermediate periods are derived. The dis- 
turbing function containing the lunar inclination factors are developed 
analytically in terms of N, the longitude of the lunar node along the ecliptic, 
rather than numerically as in reference 1. 

The results are given in tabular form making it easier to prepare 
computer programs and to select terms of given amplitude or period when 
hand calculations are desired. 

Since the development is literal, the results are applicable to three 
body systems having configurations similar to the sun, Earth, moon 
system. 
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e eccentricity of the satellite 

<»' eccentricity of the moon 

f true anomaly of satellite 

i inclination of satellite orbit 

i inclination of lunar orbit 

■t mean anomaly of satellite 

V mean anomaly of moon 

m' ratio of mass of moon to mass of earth plus moon 

n mean motion of satellite 

p summation index 

n' mean motion of moon 

q summation index 

q' summation index 

r geocentric distance of satellite 

r ' geocentric distance of moon 

T' mean longitude of lunar pericenter 
a argument of disturbing function 
a mean motion of a 

e obliquity of ecliptic 

0 argument of latitude of satellite 
0 ' argument of latitude of moon 
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LUNISOLAR PERTURBATIONS OF THE MOTION 
OF ARTIFICIAL SATELLITES 


INTRODUCTION 


For analytic purposes it is converient to start the development of the dls 
turbing function resulting from the expansion of the second order Legendre 
polynomial as given by Kozai Reference 1. The lunar disturbing function R\ 
is written in the form 


R'=n' 2 m'a 2 A *q A V c *Vq cos ^ + q ^ + v ^ 


q.q ' .v 


( 1 ) 


where 

n' mean motion of the moon 

m' ratio of mass of the moon to combined mass of moon and earth 
r radius from center of the earth to the satellite 
r* radius from center of the earth to the moon 
a, a' semi-major axis of satellite and moon respectively 
q, q v indices of summation 
c^q constants 

4) argument of latitude of the earth 
4' argument of latitude of the moon 

6 difference between the longitude of the ascending node of the earth 
and that of the moon 

i inclination of the satellite orbit to the earth's equator 
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i ' inclination of the lunar orbit to the earth's equator 
A function of i 

v q 

A', q , function of i’ 

The corresponding disturbing function due to the sun is obtained by replacing 
the prime quantities by similar quantities with the subscript © . Thus we set 

n' ~ n @ , r' = r^, a' = a Q , m' = 1 
i' ^ €, h' = 0, # = ^ 

The symbol t represents the obliquity of the ecliptic. 

The nomenclature and symbols appearing in Equation (1) are described in 
the section on symbols. 

The functions A Vq ,A' vq ,, and the constants c„ q , q are given in Table 1. The 
index q takes on the two values 0 and 2, the index q' takes the value -2, 0, 2, while 
the index v varies over the range -2, -1, 0, 1, 2. 

In certain applications, particularly when precise long range trajectories 
are required, numerical integration of the planetary equations may be useful. 

The planetary equations are therefore derived. In order to extend the accuracy 
of the computations to satellites with highly eccentric orbits, the Lagrange plane- 
tary equations are derived directly from Equation (1) without first expanding 
in powers of eccentricity. 

In order to obtain analytic formulas the disturbing function, Equation (1), 
is further expanded. After the expansion is made, the various harmonics of 
the disturbing function are shown to consist of short, intermediate, and long 
period terms. The short period terms contain the mean anomaly of the satellite. 
In the intermediate period terms the mean anomaly of the satellite is absent in 
the argument, but the mean anomaly of the perturbing body is present. For 
example, in the intermediate case for the sun, we find approximately semi- 
annual period terms, as well as multiples of this period. The corresponding 
terms of the lunar function consist of semi-monthly terms and multiples 
thereof. In the long period terms and the secular terms neither the mean 
anomalies of the satellite or the perturbing body appear. These terms have 
been treated in References 4, 5 and 6. 
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TABLE 1 



sin i 1 sin 2 i f /2 






















Since terms of intermediate period have been observed in satellite data, 
these terms are treated in this report. Analytic formulas for the perturbations 
in the motion of artificial satellites due to these terms are derived. 


THE DISTURBING FUNCTION FOR THE LAGRANGE 
PLANETARY EQUATIONS 

In order to obtain the differential coefficients needed for the planetary 
equations, it is convenient to expand Equation (1). Since q takes the two values 
0 and 2, as seen from Table 1, we obtain 


R' n' 2 m' a 2 


(jV f^rV 2Z C Vq'0 COS tvtr > 

' ' \r / (i ', v 


+ a „2 A' # C , 2 cos (2<p + q' <p‘ + vO)) 


( 2 ) 


For the convenience of the reader, the ephemeris quantities needed for 
Equation (2) are given in Appendix 1. 

The angles <j> and 0 are defined by 

~ 4 + g 

0 = h - h # 

where 

l = mean anomaly of the satellite 
g = argument of perigee of the satellite 
h = longitude of the ascending node of the satellite 
h' = longitude of the ascending node of the moon 
Some useful geometric relations, which may be derived from Figure 1, are 

cos i' = cos e cos J - sin e sin J sin N 
sin h' sin i* = sin J sin N 
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sin i' cos h 1 - sin € cos J + cos £ sin J cos N 
sin i' sin [g' - (P - N)] r. sin *. sin N 
sin i' sin J cos [g' - (P' - N)] = cos e - cos J cos i' 

v c - r 

v + g' 

■* 0' + h' 


(3) 


Equation (2) Is a form of the disturbing function useful in obtaining the 
Lagrange planetary equations. 


FURTHER EXPANSIONS OF THE DISTURBING FUNCTION 

Equation (2) while useful for forming the planetary equations is further ex- 
panded. Various characteristics of the disturbing function are then displayed 
preparing the way for analytic treatment. First, we give q’ the values -2, 0, and 
and 2 as indicated in Table 1. Then by means of Cayley’s tables, Reference 7, 
the disturbing function R’ of Equation (2) is expanded in terms of V , the mean 
anomaly of the disturbing body to obtain 


R' = n' 2 m' a 2 A *o c voo Bj 0 «• 

(i' V + 2g' + vd) 
(i' V + 24> + V6) 

+ A„ 2 A' v _ 2 c^ 22 B', -2 cos (i'V - 2g' + 24 > + v6) 

+ A V2 K 2 C v22 B 1'2 COS + 2 + 2< *> + v6 ^ 1 


+ A v0 A U C »/00 B l'0 COS + A vO A W C *20 B 'i'2 COS 

+ A ,/2 A v0 C v02 B 00 cos (20 + V0) + A v2 A v0 C v02 Bj» 0 -OS 


The coefficients B', » q ' arise from the expansion by means of Cayley's tables, 
Reference 7, and are listed in Table 2. 
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TABLE 2 


B’i'o 

i 1 

B \'-2 

i* 



B \'2 

i’ 

52/16 e’ 3 

-3 

845/48 e’ 3 

-5 

1/48 e f 3 

-1 

9/4 e’ 2 

-2 

17/2 e’ 2 

-4 

0 

0 

3/2 e’ + 27/16 e’ 3 

-1 

7/2 e* - 123/16 e’ 3 

-3 

- 1/2 e’ + 1/16 e’ 3 

1 

(1- e’ 2 )" 3 2 

0 

1 - 5/2 e’ 2 

i 

-2 

1 - 5/2 e’ 2 

2 

3/2 e» + 27/16 e’ 3 

1 

- 1/2 e» + 1/16 e’ 3 

I 

-1 

7/2 e’ - 123/16 e’ 3 

3 

9/4 e’ 2 

2 

0 

0 

17/2 e’ 2 

4 

53/16 e' 3 

3 

1/48 e» 3 

1 

845/48 e’ 3 
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Of particular interest is the quantity , representing the mean value 

of R\ The definition of <(R'^> is given by 


<R'> 


277 


I 


2 v 


R' cH 


(5) 


Since 
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<&)' 


cos (2f + /3)^ = o 2 cos /3 


where yS is an arbitrary angle, <(R')> is readily obtained. 


We find 


y ~ ^ ^ SEC + ^ | rc LP + ^ ^ INT ' 


( 6 ) 
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where R SEC represents the secular R Lp the long period and R the inter- 
mediate parts of R . 

The components of Equation (6) are, 


y ' ft sec r n<2 a2 + ~2 e ^j a oo c ooo Bqo 

£> lp = n' 2 rr # a 2 ^1 + |e 2 ^ ^T] A vO A U c v oo B oo cos ( vh “ vh ') 

+ n' 2 m'a 2 le 2 J~ > A t/2 A' v0 c„ 02 Bq 0 cos (2g + vh - vh') 

2^R' INT = n' 2 m'a 2 (l + |e 2 ^ ^] {A vO A lo c v oo B [*o cos ^h + vh') 

+ A vo A' v2 c v20 Bj, 0 cos [vh + 2 + <i' - 2) 4' - (v + 2) h'] } + 

■j n ' 2 m'a 2 e 2 J^i^ v2 A '„o c v02 B 'i'o cos (2g + vh + i'V - vh' ) + 

A v2 A v-2 c v -22 B l '-2 cos ^ 2g + - 2\' + (i' + 2)V - (v - 2) h'] 

+ A v2 A v2 c v22 B 1'2 cos [2g + I'h + 2\' + (i' - 2)V - (v + 2) h']} 


Equations (7), (8), and (9) can be used directly to find analytic perturbations 
due to the sun. 

However, as shown in Reference 1, the angle i’ appearing in the functions 
A£ q , , should be eliminated when obtaining perturbations due to the moon analytically. 
This is done in combination with angle h’ by using the first three of Equations (3). 
The quantities AJ, , are replaced by quantities A” q # p which are functions of 
the obliquity e and the inclination of the lunar plane to the ecliptic J. The argu- 
ments of the trigonometric terms will now contain N, the longitude of the lunar 
node along the ecliptic, replacing h' the longitude of the node along the earth's 
equator. 

To derive the functions A” q/p we set 


A' , cos nh' =C , 0 + C lt sin J cos N +C , 2 sin 2 J cos 2N 
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A' Vq , sin nh' = S Vq<1 sin J sin N + S Vq#2 sin 2 J sin 2N ( 10 ) 

From trigonometry we find 

A Vq' COS (Z 3 + nh ') = ^ ( C j/q* p + s# S^ q , p ) (sin J)^ P I cos (/3+p N) (11) 

p * — 2 

In Equation (11), /3 is an arbitrary angle, and 

n = - (q + v ) . 

The index p takes on the values -2, -1, 0, 1, 2. We also have the following 
relations 


for p = 0, 
for p ^ 0 


also 


c ;,<p 


- r' 

S/q'_p 


S ,'VP 


S ^-p 


O' 


r' - 1 r 

Vq' p 2 v °! p 


^vq'p 2 ^Wp 


s' = 1 fornp>0 
s' = 0 f or np = 0 

s' = -1 f or np < 0 
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We therefore find 

a ;,'p = c : q - P + s ‘ s ^p' 

C , and S tj , are given in Table 3. 

V H p vq ! ' p ® 

The quantities d Qf d, and d 2 arise from the expansion of cos' 2 i f /2 to order 
sin J. We have 


1 -2 i' 

— cos * — 
4 2 


= d 0 + dj sin J cos N + d 2 sin 2 J cos 2N 


where 


1 

2(1+ cos e cos J) 


1 sin 2 e sin 2 J 
^ (1 + cos e cos J) 2 


dj = sL n 6 

2 (1 + cos e cos J) 2 
4 _ sin 3 e 

d 2 - — 

4(1+ COS € cos J) 3 

By substituting Equation (11) into Equations (7), (8), and (9) we obtain 


A 00 A 00p C 000 B 00 (sin J)^ p ^ cos p N '2) 

P T. o 

^K'^n^m'a 2 (l +|e 2 )^] A ,o A ;o P Bi 0 (sin J) |p 1 cos (vh+P N) 

+ |n #2 m'a 2 e 2 2^A l/2 A" v0p c vQ2 Bj 0 (sin J)' p ^ cos (2g + vh + p N) (13) 


r.^SEC =n' 2 m , a 2 


1 *H 
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TABLE 3* 



The subscripts v and q' are chosen to conform with Table 1. 












_ „ > 2 „ 2 
= n * m a* 


(1 + | e2 ) Z]^oA;o p c i/oo b i-o( sin J) |P ' cos O'h + i'^pN) 


+ A^q A" c v 20 B;, 2 (sin J) 1 p 1 cos [uh + 2A' + (i' - 2)V + p NJ} + 


5 #2 #2 2 

— n ^ m a 1 * e z 
2 


J^J^vzKop c „02 B/ i' 0 (sin JV P cos (2g+vh+i +p N) + 


k vl A ;.2p C v-22 B !-2 ( sin J) P cos [2g + i^h-2V + (i' + 2)^ # + P N] 

+ A„ 2 A" v2i> c„ 22 B', 2 (sinJ) |p 1 cos [2g + vh + 2V + (i' -2)<C' + p N] > (14) 


Equations 12, 13, and 14 represent the division of the lunar disturbing func- 
tion into secular, long period, and intermediate portions. It is interesting to 
note that terms with the period of N of 18.6 years, arise from Equation (12) for 
P equal to 1 and 2. In Equations (14), i’ does not take the value zero. 

Equation (14) is used to find the analytic perturbations in the elements due 
to terms of 14 day periods or multiples of this period. 


THE LAGRANGE PLANETARY EQUATIONS 

The Lagrange planetary equations in canonical form are given by 


dL _ BR clG _ BR dH _ BR 
dt B't ’ dt ~ Bg ’ dt Bh 


d£ = BR dg = BR_ dh_ _ B_R 
dt “ BL’ dt " BG ’ dt “ BH 


The parameters a, e, and i are the semi-major axis, the eccentricity of the 
satellite and the inclination of the satellite orbit to the earth's equator, respectively. 
The Delaunay variables appearing in Equations (15) and (16) are defined by 


L = rfla , G = L /l - e 2 , H = G cos i 
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i = mean anomaly 
r - argument of the pericenter 
h - longitude of the ascending node 

We now substitute Equation (2) into Equations (15) and (16) to obtain the 
planetary equations. 

From Equation (2) let us introduce the abbreviations 


A i ~ A i/o A„ q # c Vq i 0 

A 2 = ^V2 A v q ' c vq' 2 
a i/q' = d' 0' + V 6 


We also have from Reference 8 that 





sin f 
/l - e 2 


'dl r 2 


Consequently, we find 


dL r, ,2 • 2 /a' V 

— = 2n * m — \ ^ 

Vr- j 


L\ *2>..<2* + <W ) 


sin f 


- /l - e 2 ^ A 2 sin ( 2 ^ + a v q ' H 


(18) 


dG 

dt 


= - 2 n' 2 m' a 2 'j ^ — ^ ^ | A 2 s * n ( 2 ^ + a vq > ) 


(19) 
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dH (2_* a 2 
— —— = - n m a* 

dt 


/r\ 2 /V \ s 
UJ \ r' ) 


L 


{vA, sin avq' + i/A 2 sin (2 + a ')) 


( 20 ) 


Noting that the disturbing function is written in terms of the Keplerian ele- 
ments a, e, and i it is convenient to write the planetary equations (16) in the form 


dl _ 2 /BR\ 1 - e 2 / BR \ 

dt na \Ba/ na 2 e \de/ 


dg _ cos i 

BR 

/l - | 

/ br\ 

£ ^ t na 2 e /l - e 2 

B (cos i) 

n a 2 e 

\Be / 


( 21 ) 


( 22 ) 


dh 

dt 


1 BR 

na 2 e/l (cos *>. 


(23) 


Equations (18), (19), and (20) have been derived from Equations (16) by apply- 
ing the relation 


ri 2 a 3 = /x 

da _ 2 
dL n a 


Be _ 1 - e 2 
BL _ 


na 2 e 


B e 

BG 


na 2 e 


(24) 


We now evaluate the differential coefficients appearing in Equations (18), (19), 
and (20). From Reference (9) we find 
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d 



dp 


- - cos f 


£1 
d e 



(25) 


If we make the substitutions 


B 


1 



(q' «// + v 0) 



a vJ = q' & + v 0 


we find 


) R' 

d a 


2 a 



{A„ 0 Bj + A v2 B 2 cos (20 + a^)) 


(26) 


dR' 

B e 


= -2a 2 


I) cos f £ A „ 0 B, + A„ 2 


G) 


B 2 cos (2 0 + a vq ,) 



^ — » 


sin f 7 A 2 B, sin(20+a ,) (27) 




v2 " 2 


Vq' 


BR 1 

B (cos i) 


a' 


(i )'£ 


BA 


»/0 


BA 


v2 


B (cos i) 


^ + B (cos i) “ 2 


B 2 cos (20 + a Vq# ) 


> . 


(28) 


The coefficients B A vq /B (cos i) are given in Table 1. 
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ANALYTIC PERTURBATIONS OF INTERMEDIATE PERIODS 


The method of canonical variables with the aid of a determining function, 
References 9 and 10, is used for computing perturbations of intermediate periods. 

We first write the lunar function, R*, Equation (14) in the condensed form 


R 'lNT “ n ' 2 m ' a2 A Vq A Wp B 0q B !'q' %'q ( sin J)^ ^ COS a 


(29) 


where 


a = qg + vh + + (i‘ - u) V + p N 


The complementary function R’ NT is formed by simply replacing the cosine 
terms by sine terms, to obtain 


R'.nt = n' 2 m' a 2 A„ q A*^„ B 0q B', q , c vq , q (sin J) 


sin a 


(30) 


The lunar determining function S INT is defined by 


, i V -1 r int 
>,NT ^ — 


(31) 


where 


a=qg + z'h + i'n' + p N 

The functions appearing in Equation (31) are given in Tables 1, 2, and 3. 

To obtain the solar determining function, we set p equal to zero in Equation 
(27), put A " Vq , 0 equal to A' Vq # , and make the substitution indicated below Equation 
( 1 ). 


From Reference 6, we have 

. 3nJ 2 («l+5cos 2 i) 
4 a 2 (1 - e 2 ) 2 

. -3 n J, a 2 cos i 

h = 

2 a 2 (1 - e 2 ) 2 


( 32 ) 


15 


where J 2 is the second order zonal harmonic of the earth's gravitational field, 
and a (> is the mean earth radius. 

By the theory of canonical variables we have 


b L r 


as 

a -t 


<U 


as 

dL 


as 


a k 

dG 


8H 


as 

ah 



(33) 


By substituting the expression for S given by Equation (31) into Equation (33), 
we find that the perturbations in the elements L, G, and H are 


bL = 0 

8 C. = y~" ^ INT (34) 

sh =2]^i 


The Delaunay elements L, G, and H are related to the Keplerian elements a, 
e, and i by Equation (17). 

To find the perturbations in the angular elements l , g, and h, we let 


D = n # m' A " Vcj , p B', q , c Vq#q (sin J) |P sin a. 


(35) 
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f 


Then 

U = 




( 36 ) 


where 


*B 00 

BL 


= 3(1 - e 2 ) 


L 


**02 

BL 


= 5(1 -e 2 ) 



Bg. 

?L + 


v 


Bh 

5l 


By Reference 6 


Bg _ 9 J 2 a 2 (* “ 5 cos2 O 

"5T 4 a 6 (1 - e 2 ) 2 

Bh _ 9 J 2 ag cos i 

3L 2 a 6 (1 - e 2 ) 2 


Similarly the perturbation in the agreement of perigee is 


BA 


5 g = 


_ bs _ 

B G Lmmml 


vq 


— B ftn cos i + A G 

B (cos i) 0t * v * BG 

(n /l - e 2 ) a 


3B o,' 

BG 


B 0q A V, 0 


(n /I -e 2 ) 


V 

’T 


/ 

s>. 


1 
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where 


B B, 


oo 


BG 


- 3 (1 - e 2 ) 


c')B 


02 


ac 


= -5(1- e 2 ) 


Bg 3 J 2 a e n ( 2 " 15 cos2 


BG 


2 a 2 (1 - e 2 ) 2 


<)h 15 J2 3 e n COS 1 

3G 2a 2 (1 - e 2 ) 2 


Finally the perturbation in the node of the satellite is 


'6 h = - 


as 

BH 



BA 


vq 


B cos i 
(n /l - e 2 ) a 



n /l - 




D 


(38) 


where 


Bg _ 15 J 2 a e n cos 1 

BH 2a 2 (1 - e 2 ) 2 

G Bh _ “ 3 J 2 a e n 

BH 2a 2 (1 - e 2 ) 2 
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CONCLUSION 

In this report formulas have been derived for calculating the perturbations 
on the orbits of artificial satellites due to the gravitational influences of the 
sun and the moon. Formulas 18-23 are the Lagrange planetary equations for the 
Delaunay elements which are required for predicting orbits of arbitrary 
eccentricity. 

Equations 34-38 have been derived in order to be able to calculate 
analytically solar and lunar perturbations to the orbits of satellites to within 
an estimated 0.5 meter. These equations take into account the perturbations 
of approximately 180 day periods due to the sun and the corresponding 14 day 
periods due to the moon. 
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MOON 



re 1 


APPENDIX I 


Ephemeris Quantities 


For the sun 

m 0 = .999997 

e Q = 0.01675104 - 0.00004180T 
n 0 = 0.98560027/day 

= 279f 69668 + 36000f 76892T + 0.00030T 2 
g e = 281f22083 + If 71918T + 0f00045T 2 

h © = o 

i 0 = e = 23f 452294 - 0f013013T - 0f000002T 2 

For the moon 

m' = 0.012150668 
e' = 0.054900489 

C = 270f43416 + 481267f 88314T - 0f00113T 2 
r' = 334f 32956 + 4069f03403T - 0f01033T 2 
N = 259f 18328 - 1934fl4201T 
n* = 13f 064993/day = .036291647 rev/day 
J = 5f 1453964 

where T = number of days since January 0.5, 1900, divided by 36525. 
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